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Abstract 

A model of an autonomous isothermal Brownian motor with an internal propulsion mechanism is considered. 
The motor is a Brownian particle which is semi-transparent for molecules of surrounding ideal gas. Molecular 
passage through the particle is controlled by a potential similar to that in the transition rate theory, i.e. 
characterized by two stationary states with a finite energy difference separated by a potential barrier. The 
internal potential drop maintains the diode-like asymmetry of molecular fluxes through the particle, which 
results in the particle's stationary drift. 
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1. Introduction 

Most natural (biological) molecular motors and many artificial molecular machines are powered by chem- 
ical reactions. Particular ways to convert chemical energy to mechanical work vary considerably for different 
classes of motors. A number of prototypical models were suggested which are usually rather involved and 
require application of numerical methods Q, 0, H, Hi @] ■ The purpose of this paper is to introduce a new sim- 
ple model of autonomous Brownian motor in which the complexity of involved chemical processes is hidden 
in a single parameter of an internal potential drop, not dissimilar to how an electric battery is characterized 
by its voltage. 

Our motor is essentially a Brownian particle permeable for molecules of the surrounding thermal bath, 
which can flow in and out the particle. As a physical realization of permeable Brownian particles one can 
mention for instance micro-gels Q, and lipid vesicles like liposomes When a molecule is inside the 
particle, it experiences a potential U(x) similar to that in the transition-state theory. Calculations are 
particularly simple for a linear piecewise potential U(x) depicted in Fig. 1. It is characterized, first, by the 
presence of a barrier which prevents the crossing of the particle by low-energy molecules, and, secondly, by 
a nonzero potential difference AU between the two sides from the barrier. As discussed below, the internal 
propulsion force is generated by the potential drop AU which, similar to diode devices, results in asymmetry 
of molecular fluxes through the particle and a nonzero average momentum it receives in unit time. 

In the non-operational regime, when the average molecular flux through the particle is zero, a potential 
profile like that in Fig. 1 can be formed naturally, for instance as a result of asymmetric ion distribution 
across the particle. The formation of nonzero potential difference between the two sides of asymmetric lipid 
bilayers was discussed in Q. However, it is evident that an external energy source is necessary to maintain 
an asymmetric charge distribution and internal potential drop in the operational regime of nonzero flux 
across the particle (otherwise the system would violate the second law of thermodynamics) . The function 
of a "battery" sustaining a stationary potential drop AU can be performed by particle-bound ATP- or 
light driven ion pumps, capable of moving ions against potential and/or concentration gradients (active 
transport) [jjl Hcj|. Certain proteins, which act as ion pum ps in living cells, are able to retain their function 
when transplanted into artificial membrane structures [llfl . This appears to provide a technologically feasible 
way to synthesize permeable Brownian objects with an internal potential drop. 
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Figure 1: Potential energy U(x) of molecules interacting with the permeable Brownian particle of size S. Coordinates x = 
and x = S correspond to the left and right edges of the particle. Molecules (black circles) "see" the barrier of height U a when 
coming from the left, and U a — A{7 when coming from the right. For AU > 0, the particle drifts to the right. 

The model presented here is phenomenological in the sense that specific mechanisms of active transport 
involved to sustain the potential drop AU are not identified. Of course, such an approach does not allow 
for the determination of the motor's efficiency [13], but is sufficient to evaluate transport characteristics of 
the motor as a function of AU, which is a goal of this paper. 



We consider a permeable Brownian particle of size S and mass M immersed in a thermal bath of ideal gas 
molecules of mass m <c M. The particle and molecules are constrained to move in one dimension, molecular 
distribution far from the particle is homogeneous with concentration n, and the velocity distribution of 
molecules is Maxwellian 



where vt — 1 / \fmfi and j3 — I/HbT. Molecules interact with the particle with a piecewise linear potential 

U{x) depicted in Fig. f. The direct interaction of molecules to each other is assumed to be negligible. 
Suppose the particle's interior corresponds to the interval < x < S, with x = So being the position of the 
barrier's maximum. When inside the particle, a molecule experiences a constant force / = \dU/dx\ directed 
to the left on the left slope < x < So, and to the right on the right slope So < x < S. The right arm of 
the potential is shorter than the left one by AS, which results in the difference of baseline potential values 
(potential drop) AU = fAS. 

The idealized model potential like that in Fig.f can be formed by two electric double layers with opposite 
polarities, infinitely extended in the y and z directions. In a realistic case of the particle of finite dimensions, 
one has to take into account a nonzero field outside the particle and, in general, a discontinuity of the 
potential at the particle boundary [l3j], but these would not change the problem qualitatively. 

It is intuitively clear that for given geometry the average force exerted on the particle by molecules 
coming from the left exceeds the force from molecules coming from the right, and thus the particle would 
drift to the right. Our aim is to find the stationary drift velocity of the particle (V). (A more simple 
model with a monotonic step-like potential U{x) shows a nonzero drift, but no friction, and therefore, in 
our opinion, is of less interest.) 

In lowest order in the mass ratio parameter m/M , which is sufficient for our purposes, one can justify 
the standard Langevin equation for the particle's velocity V , 



2. Model 




(1) 




M 



dV(t) 
dt 



1 V{t)+F(t), 



(2) 



where however the fluctuating force F(t) is not zero-centered, so that the average stationary drift velocity is 



Even for the simple model outlined above, the analytic evaluation of (F) and 7 is rather involved except 
when the potential drop is small and the barrier is high, 



/3AC/«1, pU a >l. (4) 

The first condition allows one to evaluate the average fluctuating force in a linear approximation in AU, in 
which case the drift velocity can be written as 

(V) = aAU, (5) 

where the parameter a can be called the mobility. The second condition, (3 U a > 1, implies that the friction 
coefficient 7 can be calculated in the limit of an infinitely high barrier, when the fluctuating force F(t) is 
zero-centered, and the standard fluctuation-dissipation relation holds 

/>oo 

7 = /?/ dt(F(0)F(t)). (6) 
Jo 

From a practical point of view the conditions Q are of less interest since the corresponding mobility a 
is very small. Yet in what follows we evaluate a assuming for the sake of simplicity that conditions (U) are 
satisfied. One naturally expects that in this case the temperature dependence is given by an Arrhenius law 
er = Co exp(— j3U a ), so that our major interest will be in the pre-exponential factor ctq- 



3. Impermeable particle 

In this section, wc evaluate the correlation function of the Langcvin fluctuating force F(t) and, using the 
fluctuation-dissipation relation ((SJ), the friction coefficient 7, for the case of a high barrier or low temperature, 
when the molecular barrier crossing is negligible. Clearly, in this case of an impermeable particle the potential 
drop does not play any physical role and can be set to zero, AU — 0. Recall that we need 7 in order to find 
the drift velocity of the motor with Eq.©. As will be found below, for the given model 

[2 

j = 4y — nmvT, (7) 

where n is the concentration of bath molecules. Equipped with this result, the reader not interested in 
derivation details may skip the rest of this section. Yet, since our evaluation is microscopic and exact, it is 
also a problem of interest of its own. It is also not entirely distracting since in the next section a similar 
method will be applied to evaluate (F) for a permeable particle. 

It is convenient to write the fluctuating force as the sum of contributions from molecules coming from the 
left and right, F = Fl + Fr. For the given case F is zero centered, (F) = 0, and the statistical properties 
of Fl and Fr are the same. Then it is sufficient to consider only one, say the "left" , contribution, since the 
correlation function of the total force can be written as 



(F(0)F(t))=2((F L (0)F L (t))), (8) 

where double brackets denote the cumulant ((A B)) = (AB) — (A) (B) . 

As known from the microscopic theory of Brownian motion [16l |. to lowest order in the small mass 
ratio parameter m/M the fluctuation force F(t) can be evaluated neglecting the particle's motion, i.e. as a 
pressure on an infinitely heavy particle. Suppose the particle's left edge is fixed at x = 0, so that i-molecules 
(coming from the left) experience the potential U{x) which equals fx for x > and zero for x < 0, see 
Fig. 2. Each L-molecule inside the particle (x > 0) gives the same contribution / to the total force on the 
particle, which therefore can be written as 

/oo y* 00 

dv / dxN(x,v,t), (9) 
-00 Jo 
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Figure 2: Left figure: a snapshot at time t of three L-molecules with the same initial velocity v which are inside the particle 
(x > 0). Arrows show velocity directions. Molecule 1 has just entered the particle, molecule 2 has reached a maximal position, 
molecule 3 is about to leave the particle's interior. Right figure: a snapshot of the same molecules at an earlier time t — r r . All 
three molecules are outside the particle and occupy the interval — vr r < x < 0. 



where 

N(x,v,t) =^28(x-Xi)5(v-Vi) (10) 

i 

is the microscopic phase density of L molecules. 

Let us define the residence time T r (v) as the time which an L-molecule with initial velocity v spends 
inside the particle, i.e. in the region x > 0. Suppose a molecule enters the particle at t — with the initial 
velocity v. Inside the particle the velocity changes with time as v(t) = v — (/ /m)t. At the moment when 
the molecule leaves the particle, its velocity is —v. Then the equation v — (//m)r r = — v gives the residence 
time 

. „ 2mv 

Tr(v) = — • (11) 

For brevity, we shall often suppress the argument v of the function T r (v). 

In order to evaluate the correlation function of the force Fi,(t), it is convenient to express it in terms of 
phase density at an earlier time t — T r . All molecules, which at time t are inside the particle x > 0, at time 
t — 7> were outside in the region — vr r < x < 0, see Fig. 2. Then, instead of one can write 

r0 



F L(t) = f dv dxN(x,v,t-T r ). (12) 

JO J-VT r 



The advantage of this form is that it refers to molecules in the potential free region x < where the velocities 
of molecules are constant, so that the following relation holds 

N(x, v, t) = N(x - vt, v, 0) = N {x~ vt, v), (13) 

where the notation Nq(x, v) is introduced for the density at time t = 0. This relation allows one to reduce 
temporal correlations to single-time correlation in phase space. For an ideal gas the latter takes the form 

({N (x,v)N (x',v')}} =nf M (v)5(v-v')6(x-x'), (14) 

where Jm{v) is the Maxwellian distribution ([T]), and n is the density of the gas. 

With the above relations, the evaluation of the correlation function for the force F^t) can be done as 
follows. From (fl^)) . one gets 

/>oo />oo />0 p0 

((F L (t)F L (0))) =f \ dv dv' dx dx' ({N(x,v,t- Tr )N(x', «',-<))), (15) 



JO JO J—vr r J —v' t' t 

where t' t — T r (v'). Using (|13|) . we can write the cumulant of the phase density which appears here as 

((N {x -tv + T r v,v)N (x' + T' r v',v'))). (16) 
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From ([H )l -(|16 p , after integration over v r , one obtains 

POO 

((F L (t)F L (0))) =nf 2 dvf M (v) I dx I dx' 6(x - 1 v - x'). (17) 



The integral of the delta function over x' equals one when 

- vr r < x -tv < 0, (18) 

and zero otherwise. Since Eq. pTl) involves only negative x, the right side of the condition (1T51) is satisfied 
(we assume t > 0), while the left side requires x > —v(r r — t). Then Eq. dl~7l) takes the form 



((F L (t)F L (0))) =nf 2 / dvf M (v) dx0(x + v{T r ~t))=nf 2 dv f M (v)v {r r - t) 9(r r - 1), (19) 

JO J-vt t JO 

where 0(x) is the Heaviside step-function. Recalling expression (jlip for the residence time, one obtains 



{{^.|/).F i (())}}=2,//„/ / ( /r/,;(r) (20) 



/t/2m 

which eventually gives 

(F(t)F(O)) = 2((F L (t)F L (0))) = 2nm/4 erfc , (21) 

where erfc (a;) = 1 — erf (a;) is the complementary error function, and the characteristic time of correlation is 
to = \/2T r (vT) — 2y/2m Vt/ f- Then for the friction coefficient 7 we obtain the result (JT]), 

7 = /3^ dt{F(t)F(0)) =4^nmv T . (22) 

It is interesting to note that this expression, which we found here for a linear potential U (x) (constant 
repulsive force), coincides with 7 for a similar model where molecules interact with the particle with a 
parabolic potential (linear repulsive force) ■ The same result for 7 can be obtained from the model with 
instantaneous (and therefore binary) particle- molecule collisions 15( . For these models correlation functions 
(F(0)F(t)} are different, but the integral (j2"2l is the same. 



4. Permeable particle 

In this section we assume that the potential barrier height is large but finite, /3U a > 1, so that hot 
molecules of the thermal bath can pass through the particle. For a non-zero internal potential drop, AU 7^ 0, 
the fluctuating force F in the Langevin equation ([2]) is not zero-centered, which causes the particle's drift. 
Our goal is to evaluate (F) and the particle's drift velocity (V) = (F)/j (with 7 found in the previous 
section) within a linear approximation, treating AU as a small perturbation. 

As in the previous section, we assume that the particle occupies the interval < x < S, and we write 
the force as the sum of contributions from molecules coming from the left and from the right, F = Fl+ Fr. 
Consider first F^. It is convenient to analyze it as the sum of two parts, 

F l =f£ ] +F {2) . (23) 

The first part i^ 1 comes from "slow" molecules with initial velocities (before entering the particle) less 
than the activation velocity 



(24) 




Figure 3: Left figure: a snapshot of four fast L- molecules with the same initial velocity v > v a . Only molecules located in 
the interval < x < AS (molecules 1 and 2) give nonzero contribution to the average force (Fl) exerted on the particle by 
L-molecules. Contributions from molecules with AS < x < So (molecule 3) and So < x < S (molecule 4) cancel out. Right 
figure: a snapshot of the same molecules at an earlier time t — At, when molecules 1 and 2 were are outside the particle. 

(2) 

which is the minimal speed for L-molecules to cross the barrier. The contribution ' corresponds to "fast" 
molecules with initial velocities v > v a . 

For molecules with v < v a the barrier is impenetrable, and in order to evaluate one can borrow the 
results of the previous section. Namely, as illustrated in Fig. 2, one can write i<2 at a given time t in terms 
of the phase density of molecules at an earlier time t — r r as follows 



Here the residence time r r (v) is given by (fTI]) . and / = \dU/dx\. For the potential-free region (N(x,v,t)) = 
nfjuiv), then from (|2"5"j) and ([TTj) one obtains 



In the limit of impermeable particle, v a — > oo, the above equation gives the result of the elementary kinetic 
theory for the pressure force (F) — mn(v 2 ) — n//3. 

Consider now the second part F L in Ea. ((23|) . which comes from fast molecules with initial velocities 
v > v a . Such molecules cross the barrier and can occupy both sides of the potential hill. Clearly, the average 
forces exerted by fast L-molecules on the left and right slopes cancel out, except for the molecules near the 
base of the left slope with coordinates < x < AS, see Fig. 3. Let Ar be the time for an L-molecule with 
initial velocity v to climb to the altitude AU, 




(25) 




(26) 




(27) 



Similar to (|25p , one can write 




(28) 



which gives 




(29) 



The same method can be applied to evaluate two parts of the force 



Fn - L (1) + F 



R ' 



(30) 



G 



exerted on the particle by i?-molecules. The only difference is that for molecules coming from the right the 
potential barrier is lower, U a — At/, and so is the activation velocity v a = [2(U a — AC/)/m] 1 / 2 . In the linear 
approximation, 

~ a a V a&U 

v a ^v a -Av a , Av a = , (31) 

where v a is the activation velocity for the L-molecules given by Ea. (|24j) . Using the same argument as for 
L-molecules, the contribution from slow i?-molecules with |u| < v a , which cannot cross the barrier, can be 
written as 

(F^ ) )=-2mn dvf M {v)v 2 . (32) 
Jo 

In the linear approximation this takes the form 

(F$p) = -2 m n / dv v 2 + 2mn ,fM{v a )v 2 a Av a , 



or 

(F^) = ~(F^) + 2nf M (v a )v a AU. (33) 

The contribution from fast i?-molecules with \v\ > v a is similar to Eq. (|2T)]) for fast L-molecules. In fact, in 
the linear approximation the two expressions coincide, 

/>OC POO 

(F {2) ) = nAul dvf M (v)^nAU dv f M (v) - (F (2) ) (34) 
Combining the above results, the total average fluctuating force can be written as 

(F) = \2nf M {v a )v a + 2nj dv AU. (35) 

Explicitly this reads 

(F)=n4>(^J AU, (36) 

where the dimensionless function <p(x) is 

#*0 = ^exp(-^) + erfc(^). (37) 



Here the first and second terms originate from slow and fast molecules, respectively. For large x (high 
barrier/low temperature) the second term is smaller than the first one by the factor of order 1/x 2 and can 
be neglected. Then, since v a /vT = (2/3C/Q) 1 / 2 , one gets 

2 



(F) = —n ^fWl e- 0u ° AU. (38) 
V 71 " 

This and the expression Q for 7 give for the stationary drift velocity (V) = (F) of the motor the following 
result 

1 n lUa 



(V)=*AU, a= — ^ e -^. (39) 

As expected, the temperature dependence of the mobility a is Arrhenius-like. The square-root dependence 
of prc-cxponcntial factor on the barrier's height U a is perhaps a less obvious result. The mobility a depends 
neither on the mass of the motor nor on the concentration of molecules of the thermal bath. 
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Figure 4: Numerical simulation of the drift characteristics of the motor with /3U a = 5, /3A{/ = 0.05, and mass ratio 
m/M = 0.1. Upper figure: average displacement of the motor as a function of time. Lower figure: average drift velocity. Solid 
lines correspond to numerical experiment data, averaged over 1.8- 10 7 trajectories. Dashed lines represent theoretical predictions 
given by Eg . 14011 . Units of velocity, time, and displacement are, respectively, no = vt, to = T r (vT)/2, and xo = voto- 

5. Discussion 

The autonomous mesoscopic motor described in this paper is based on a permeable Brownian particle 
which propels due to asymmetric transfer characteristics, and in this sense may be called a Brownian diode. 
The main result (|39p for the motor's mobility a is obtained under asymptotic conditions of a high barrier 
and small potential drop To verify this result, we performed a one-dimensional numerical simulation 
similar to standard molecular dynamics with the velocity Verlet algorithm. Instead of periodic boundary 
conditions that poorly fit the problem, a very large thermal bath with a Maxwcllian velocity distribution 
([T]) and a constant density of molecules was generated in the beginning of each simulation run. The size 
of the bath was of order cvxts, where t s is the duration of one simulation run, and c is a numerical 
parameter of order 10 2 . As far as the Brownian particle dynamics is concerned, the difference between a 
bath of such size and the infinite one is that the particle immersed in the former may feel a deficit of hot 
molecules with velocities higher than cvt- For c > 10 2 , it was found empirically that results do not depend 
on c, and thus the difference with the infinite bath is expected to be negligible. The assumption of non- 
interacting bath molecules strongly reduces the amount of computations. It also allows some optimization 
by eliminating molecules with initial velocity directions away from the particle and also those which are so 
slow and generated so far from the particle that cannot come to its proximity during the time of simulation. 
The ensemble average of the particle's coordinate and velocity as functions of time is taken by performing 
multiple runs, resetting the bath molecules with new initial conditions and averaging over the simulation 
runs. 

Simulation shows a good agreement with the prediction for small AU and when f3U a > 5, see Fig. 4. 
The agreement gets better when one takes into account the contribution of fast molecules given by the 
second term in Eg. ((37)) . Since erfc(x) rj e~ x / \fiix for large x, the corrected expression for the mobility 
takes the form 



For lower barrier heights, j3U a < 5, the simulation shows that Eqs. (|3U)) and (|40[) both noticeably un- 
derestimate the mobility, see Fig. 5. The discrepancy comes from the approximation of an infinitely high 
barrier, which we have used to calculate the friction coefficient 7, Eq. ([22)) . Clearly, this overestimates 7 and 
underestimates the drift velocity (V) = (F)/j. An analytical evaluation of 7 for the case of a finite barrier 
is feasible but somewhat cumbersome. One can instead find 7 from a numerical experiment measuring 




(40) 
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Figure 5: Same as Fig. 4, but for a lower potential barrier (3U a 
theoretical results becomes more significant. 



4. For a lower barrier the disagreement of experimental and 



velocity relaxation (V(t)) = Vq cx-p(-jt). When one uses 7 found in this way, and (F) given by a good 
agreement with the numerical experiment may be restored. (This method is not very effective for the mass 
ratio m/M > 0.1, when relaxation of the Brownian particle's velocity is distinctly non-exponential.) 

In our calculations, we assumed that the two slopes of the potential hill are equally steep. No qualitatively 
new feature appears if they are not (provided the potential is continuous). As shown above, the stationary 
mobility does not depend on U'(x). Respectively, for an asymmetric barrier with no potential drop, the 
drift may be only transient and vanishes in the long-time limit. This is confirmed by numerical simulation 
for regimes within and beyond the linear response. At this point, however, it perhaps should be recalled 
that in this study we ignored the particle's internal degrees of freedom. Taking them into account, one may 
expect that the model would reveal new features characteristic for granular Brownian motors, interacting 
with bath molecules via inelastic collisions [l7j |. 
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